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Abstract 

We investigate two important questions about the use of the nonextensive thermo- 
statistics (NETS) formahsm in the context of nonhnear galaxy clustering in the 
' Universe. Firstly, we define a quantitative criterion for justifying nonextensivity 

■ at different physical scales. Then, we discuss the physics behind the ansatz of the 

! entropic parameter q{r). Our results suggest the approximate range where nonex- 

' tensivity can be justified and, hence, give some support to the applicability of NETS 

^ ■ to the study of large scale structures. 
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1 Introduction 



The evolution of large scale structures in the Universe is one of the most 
important questions of modern cosmology. A considerable amount of work 
has been done on numerical and analytical approaches to characterise the 
clustering of matter at large scales. One of the difficulties inherent to this 
subject is that, in order to make progress in the understanding of general 
models, it is necessary to define methods for structure quantification. 

Usually, quantitative studies of large scale structures are based on the two- 
point correlation function ^(r) for the galaxy distribution. Estimates indicate 
this function is well approximated by ^(r) = (r/ro)"''' (with 7 ^ 1-8), where 
tq is the correlation length, the scale marking the transition between lin- 
ear and non-linear regimes. The usual range for 7, found in the literature, 
is 1.5 < 7 < 1.97 (see, e.g., [1] [2]), with tq = 5h~^ Mpc for galaxies and 
ro = 25h~^ Mpc for cluster (see, e.g., [3]). At larger scales, the structure of the 
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Universe presents patterns like walls and filaments (with dimensions ~ 1501i~^ 
Mpc) seeemgly reaching homogeneization at Cosmic Microwave Background 
(CMB) scales lOOOh"^ Mpc). Some authors have had some success describ- 
ing the clustering properties of visible matter over this wide range of scales in 
terms of a multifractal phenomenon associated with density thresholds applied 
to multifractal sets (e.g. [4] [5]). However, the relative success of the multi- 
fractal approach does not imply a better understanding of the physics behind 
this framework. Actually, it is not simple to find a dynamical connection be- 
tween fractal sets and galaxy clustering. Recently, Ramos et al. [6] (hereafter 
RWRR) shed some light on this question by putting forward a model based 
on the generalized thermodynamics (NETS) formalism [7]. They show that 
applying the idea of nonextensivity, intrinsic to NETS, it is possible to derive 
an expression for the correlation function 

1 + e(r) = ^D2r(^-^\ (1) 

using a scale-dependent correlation dimension 

logP + oO^ (2) 
log 2 



where the entropic parameter q{r) is given by the following ansatz 



with a and /3 being free parameters of the model. This approach shows a 
smooth transition from a clustered Universe to large-scale homogeneity, with 
D2 = 3. However, RWRR do not discuss two important questions concerning 
the conceptual basis of the model: the necessity of defining a criterion which 
allows us to assume nonextensivity in the context of galaxy clustering and the 
physics behind the ansatz for q{r). In the next two sections these questions 
are further developed and some conclusions are drawn at the end. 



2 Nonextensivity and gravitational clustering 



The physical motivation of the approach adopted by RWRR is built upon the 

fact that components of gravitating systems tend to evolve spontaneously into 
increasingly complex structures due to the long range nature of the gravita- 
tional interaction. The NETS theory generalises the Boltzmann-Gibbs statis- 
tical mechanics and can be applied to systems dominated by the long-range 
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nature of gravity. However, the application of NETS to an ensemble of co- 
moving cells containing gravitating particles depends on the behaviour of the 
average correlation energy inside a spherical cell of volume V = iirR^/S with 
increasing scales R. For instance, the grand canonical ensemble of cells that 
are larger than the correlation length are approximately extensive, since the 
universal expansion effectively limits the thermodynamic effects of gravity to 
roughly the correlation length scale (e.g. [3] [8]). Also, for an individual cell 
whose size is larger than the correlation length, extensivity is possibly a good 
approximation because the correlation energy between two members of the 
ensemble is negligible compared to the internal correlation energy: 

Ucorr < t/l + t/2 ^ Um ^ Ui + U2. (4) 



Due to these reasons, the application of NETS to galaxy clustering is not 
straightforward and demands better criteria to properly describe the prob- 
lem using the non-extensive formalism. For real nonextensive systems, the 
correlation energy between two cells should be as important as the internal 
correlation energy, such that 

Utot = C/l + C/2 + Ucorr- (5) 



Following Sheth & Saslaw [9] , we compute the average gravitational correlation 
energy within cells of volumes V and 2V for increasing R: 



CirUTTfir'^ dr, (6) 





where m is the mass and n is the average number density of particles in 
a cell of size V. The extensivity approximation requires that W2V ~ 2VFv, 
approximately verified for the power law correlation, ^(r) = ^o'""''', whenever 
7 > 1 (see [9]). In this case, |1^2v/2VI/V| = 2*^^^"'''/'\ which means that, using 
7 = 1.77, we find a ~5% deviation from the strict extensivity condition at all 
scales. 

In the NETS context, we substitute (1), (2) and (3) into (6), numerically in- 
tegrate for V and 2V (the spherical cell of volume 2V has radius 2^^^R) and 
calculate the correlation ratio \ W2v/2Wv\ for each upper limit R. The results 
are presented in Figure 1, where we plot the three cases investigated by RWRR 
in comparison to the line defined by the power law correlation (for 7 = 1.77). 
Note that, for the three NETS models, we see different levels of deviations 
from the extensivity approximation. In particular, NETS- 3 model presents 
deviations from extensivity of ~20% even at very large scales. However, for 
large enough scales, the correlation ratio decreases and the extensivity approx- 
imation is recovered, for models NETS_1 and NETS_2. In order to properly 
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invoke noncxtcnsivity, we propose a lower limit for the ratio \W2v/^Wv\ as 
10%. This is twice the energy ratio for the power law correlation case. 
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Fig. 1. Behaviour of the correlation energy ratio with the scale R. 

One can further quantify that as follows. Consider two elliptical galaxies of 
radii Vg separated by the distance tq. A family of analytic models for spheroidal 
stellar systems is defined by the density distribution: 



V 



47r r3-'?(l + ry+v 



(7) 



for < < 3 [10]. Choosing units in which the total mass and the gravitational 
constant are both unity, we find the self-gravitational energy of each galaxy 
as 



2 \2r]-l^ 



(8) 



and the potential at a distance Tq is 
1 



77-1 



(1 + ro)"-! 



(9) 



with the mass interior to radius rg being 



MM 



{l + rg)v- 



(10) 
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The gravitational energy of the interacting system formed by the two galaxis 
is given by 

Uint^'2%{ro)Mr,{rg). (11) 



Assuming we have only the two galaxies in the volume defined by AttTq/S, we 



should have Ur, 



Uint- Then Utot = 2^7„ + Uint and, consequently, 



tot 



2Ur, 



2r]-l 



.(l + ^o)"- 



[12) 



In Figure 2 we present the behaviour of Utot/'^U^j as a function of Tq (for 
ro < 5h~^ Mpc and taking rg — 30h~^ kpc). Note that only at very small scales 
the ratio significantly increases. Actually, at the galaxy correlation length, the 
deviation from extensivity is about 1%. As a comparison, NETS models give ~ 
13%, 17% and 28%, for the cases 1, 2 and 3, respectively, at the same scale, in 
the more general situation. Thus, a NETS model with at least 10% of deviation 
from the strict extensivity condition corresponds to 10 times the deviation in 
the case of a cell with only two typical elliptical galaxies interacting within 
the galaxy correlation length (ro = 5h~^). 

Using this criterium, NETS_1 and NETS_2 models are well justified for r < 10 
h~^ Mpc and r < 30 h~^ Mpc, respectively, while NETS_3 model is valid at 
all scales. This result reinforces the idea of applying the NETS formalism to 
galaxy clustering phenomena, but clearly indicating the approximate range it 
will be used, given a specific choice of the free parameters a and /3. Hence, the 
level of nonextensivity strongly depends on the the ansatz (3), which reinforces 
the need of a deeper discussion of the meaning of q{r). 



3 The physics behind q(r) 



Structure formation models try to compute the evolution of cosmic structure 
from the very early Universe to the present day. Usually, they are stochastic, in 
the sense that random initial conditions are used, with well-specified statistical 
properties, and their late evolution may be highly non-linear [11]. Any attempt 
to compute g(r) from first principles must take into account this complex 
scenario. 

Considering the relevant role played by turbulence in the dynamics of struc- 
ture formation (e.g. [12] [13]), a natural approach for deriving a meaningful 
expression for g(r) is to assume a. a la Kolmogorov cascade phenomenology. 
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Fig. 2. Correlation energy ratio for r/ models. For comparison, the Hernquist's model 
is computed for ij = 2. 

in which the energy supplied at large scales by the physical mechanisms act- 
ing over structure formation flows down until being flnally dissipated on the 
smallest scales by viscous processes [14] . The key ingredient in this "top-down" 
cascade scenario is the presumption of the existence, within a certain range 
of scales, of a scaling (w") ~ r^" of the moments of the peculiar velocity 
differences Vr{x) = v{x + r) — v{x). This scale-invariance can be rigorously 
deduced, under assumptions such as local isotropy, from the dynamical equa- 
tions governing a turbulent fluid [15]. In the present cosmological context, 
scale-invariance is well supported by observations [16] and is at the heart of 
the fractal description of galaxy clustering [4] . If we now assume that the prob- 
ability density function (PDF) of peculiar velocity differences is described by 
the Tsallis canonical distribution and that at sufficiently large scales turbulent 
fluctuations are normally distributed, then an analytical expression for q{r) 
can be easily derived [17], in the form of g ~ (15 — 21r") / (9 — 15r"), where 
a = (4 — 2(2 and (p are the structure function exponents. 

This model was successfully applied to hydrodynamics turbulence ([18], [19]), 
where the entropic parameter represents a direct measure of intermittency. It 
gives similar results to those obtained with equation (3), albeit only within a 
restricted range of scales. Note that, as r — 0, the entropic parameter tends 
to a flnite value, which is physically expected for any model of q{r) [20]. But in 
this case, equation (3) would not provide a smooth transition from small-scale 
fractality to large-scale homogeneity, which represents a major drawback in 
the present cosmological context. 



The problem with such a cascade phenomenology is its description of structure 
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formation from larger to smaller scales. To support observations, a different 
and more compelling model may be obtained from a "bottom-up" fractal cas- 
cading scenario. In this cold dark matter (CDM) dominated Universe scenario, 
large scale structures are formed by gravitational clustering of smaller clumps 
of matter [20]. This process eventually leads to strong density fluctuations 
with self similar density distribution and a stationary fractal dimension, while 
enhancing long-range correlations and the corresponding high energy tails in 
the peculiar velocity differences PDFs. 

In this context, averaging a Gaussian conditional velocity distribution over all 
possible spatio-temporal energy dissipation rate fluctuations and at appropri- 
ate scales leads to a Tsallis peculiar velocity differences PDF and to a closed- 
form expression for the entropic parameter [21]. This expression depends on 
the (discrete) number of degrees of freedom relevant to represent the local fluc- 
tuations. It can be further generalized, taking the form g ~ (r -|- 3) / {r + 1). 
Considering the lack of observational data for an unambiguous determination 
of q{r), we foresse a complementary (and computationally expensive) approach 
to gain a deeper insight on the role of the entropic parameter. We propose to 
compute estimates of q within a large range of scales directly from A'^-body 
simulations, using a A — CDM cosmological model of gravitational clustering. 



4 Summary and conclusions 

The basis for clustering statistics in cosmology is the study of galaxy distribu- 
tion in the Universe. The analysis of an increasingly amount of astronomical 
data can provide the correct framework to explain the formation and evolu- 
tion of large-scale structure in the Universe. In particular, correlation functions 
works suggest the possibility of describing the galaxy clustering properties in 
the context of the multifractal approach. Recently, RWRR presented a bet- 
ter physical interpretation to this approach by deriving multifractality from 
the NETS formalism. Now, extending their work, we define a quantitative cri- 
terium to use the NETS formalism based on deviations from strict extensivity. 
Assuming it should reach at least 10 times the expected deviation level within 
a cell containing only two typical galaxies, we obtain the approximate validity 
domain for the NETS models: r < 10 Mpc for NETS_1 model, r < 30 Mpc 
for NETS_2, and at any scale for NETS_3. The choice of q{r) has a physical 
meaning in the context of hydrodynamics turbulence and can be understood 
as a "bottom-up" fractal cascading, in conceptual agreement with CDM struc- 
ture formation models. Although the scenario we propose is not complete, we 
should keep in mind that there is not a unified framework to explain galaxy 
clustering in the Universe. Further numerical developments of NETS models 
will be presented in a forthcoming paper [22]. 
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